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We present results of a high statistics study of the quenched spectrum using Wilson 
fermions at /3 = 6.0 on 32^ x 64 lattices. We calculate the masses of mesons and baryons 
composed of both degenerate and non-degenerate quarks. Using non-degenerate quark 
combinations allows us to study baryon mass splittings in detail. We find significant 
deviations from the lowest order chiral expansion, deviations that are consistent with 
the expectations of quenched chiral perturbation theory. We find that there is a ~ 20% 
systematic error in the extracted value of mg, depending on the meson mass ratio used 
to set its value. Using the largest estimate of rus we find that the extrapolated octet 
mass-splittings are in agreement with the experimental values, as is Ma — Mat, while the 
decuplet splittings are 30% smaller than experiment. Combining our results with data 
from the GFll collaboration we find considerable ambiguity in the extrapolation to the 
continuum limit. Our preferred values are M^/Mp = 1.38(7) and M/^jMp = 1.73(10), 
suggesting that the quenched approximation is good to only ~ 10 — 15%. We also analyze 
the 0{ma) discretization errors in heavy quark masses. 
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1. Introduction 

Precise measurements of tlie liadron spectrum using lattice QCD are crucial both to 
validate QCD as the correct theory of strong interactions and to establish the reliability 
of numerical simulations for extracting weak matrix elements. Current lattice calculations 
suffer from a variety of systematic errors, most notably those due to quenching, discretiza- 
tion and the need to extrapolate to light quark masses. In this work we present a detailed 
study of these systematics for the hadron spectrum. 

Such a study requires small statistical errors. We have reduced these by using a 
moderately large ensemble, 170 lattices, and working on a large lattice, 32^ x 64 at = 6.0 
in the quenched approximation. We use unimproved Wilson fermions. Preliminary results 
from a subset of 100 lattices were presented at the LATTICE94 meeting [0. Our lattices 
are large enough that we expect finite size effects to be small. The major technical features 
of our work are (a) using two kinds of sources that yield correlators that converge to their 
asymptotic values from opposite directions, so as to improve the reliability of the masses 
extracted; (b) calculating hadron masses using non-degenerate light quarks, which allows 
us to study the quark-mass dependence in detail; and (c) using ratios of correlators to 
obtain accurate estimates of mass differences. We find that terms of higher order than 
linear in the quark mass are very significant, and that their inclusion is essential for the 
extrapolation to the physical light quark masses, particularly for mass splittings amongst 
spin-1/2 baryons. Higher order terms are also important for the vector mesons and the 
spin-3/2 baryons. 

The outline of the paper is as follows. In the following section we explain how we 
generate lattices and calculate quark propagators. After a brief discussion of fitting, and a 
summary of the expected chiral behavior of hadron masses in the quenched approximation, 
we present our results for mesons and baryon masses. We then extrapolate mass ratios 
to the continuum limit by combining our results with those of the GFll collaboration 0. 
There turns out to be considerable ambiguity in this extrapolation. Our preferred values 
for the extrapolated ratios are /Mp = 1.38(7) and Mi^jMp = 1.73(10). This suggests 
that the quenched approximation is good to ~ 10 — 15%, less accurate than suggested in 
Ref. 1^. We close with some conclusions and suggestions for further work. 

We use the following conventions throughout the paper. Hadron masses are denoted 
by upper case M, while for quark masses we use lower case m. All masses are in lattice 
units unless explicitly expressed in MeV (or GeV). 



2. Details of simulations 

Gauge configurations are generated using a combination of over-relaxed (OR), pseudo- 
heatbath (PHB) and Metropolis algorithms. Typically, 5 OR sweeps are followed by a PHB 
update, with the latter consisting of three hits, one in each of the SU(2) subgroups. In some 
cases the PHB update is replaced by a 20-hit Metropolis sweep. Two independent streams 
were generated, each starting from a lattice consisting of two independently thermalized 
32^ lattices joined together. A further 2000 x (50R + IPHB) thermalization sweeps were 
then performed. Thereafter we analyze lattices separated by 400 x {50R+1PHB) sweeps. 
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We calculate quark propagators using the simple Wilson action, with periodic bound- 
ary conditions in all 4 directions. We use two kinds of extended sources - Wuppertal and 
Wall - at each of the five values of quark mass given by k = 0.135 (C), 0.153 {S), 0.155 
(Ui), 0.1558 (f/2), and 0.1563 {U3). These quarks correspond to pseudoscalar mesons of 
mass 2835, 983, 690, 545 and 431 MeV respectively, using 1/a = 2.33 GeV for the lattice 
scale. We use the three light quarks to extrapolate the data to the physical isospin sym- 
metric light quark mass fn = {rriu + md)/^, while the C and 5" k values are selected to 
be close to the physical charm and strange quark masses. The physical value of strange 
quark, in fact, lies between S and Ui and we use these two points to interpolate to it. In 
most cases we find that the extrapolation to m can be done using the six combinations of 
light quarks UiUi, U1U2, UiUs, U2U2-, U2U3, U3U3. For brevity we use {UiUj} to refer 
to this set of masses. 

We analyze three types of hadron correlators distinguished by the sources and sinks 
used to generate quark propagators. These are wall source and point sink (WL), Wuppertal 
source and point sink (SL), and Wuppertal source and sink (SS). The notation and details 
of the implementation of the Wuppertal source are as in our previous work 0. The 
smearing parameter is set to i^kg = 0.181, corresponding to a smearing size of ~ 3. 
(In this was mistakenly written as O ~ 3.) The generation of the Wuppertal sources is 
a negligible overhead on the inversion. 

The Dirac equation is solved using the over-relaxed preconditioned (fourth order poly- 
nomial) minimal residue (MR4) algorithm described in Ref. [Q. The convergence rate is 
quite insensitive to the over-relaxation parameter u as long as it is in the range 1.2 — 1.35; 
we use uj = 1.3. We set the convergence criteria (for all values of k) to < 10~^^, 

where r is the remainder and x is the solution. This tolerance is as small as we can de- 
mand, as we use IEEE single precision arithmetic. To ensure convergence we run the MR4 
inverter up to three times, refreshing the starting remainder each time, and then switch 
to conjugate gradient, and force it to run at least one cycle. To date, we have observed no 
failures of MR4. We note that the simple MR algorithms are much less sensitive to our 
use of 32-bit arithmetic than the CG-based algorithms, whose convergence rate depends 
on maintaining orthogonality of a sequence of vectors. Indeed, we have found that MR4 
is the most efficient and stable of the algorithms (MR, MR2, MR4, CO, BiCG5 [g and 
BiCGStab [0) we have implemented on the CM5. 

A technical detail of our MR4 algorithm that makes it suitable for 32 bit precision is 
as follows. The n*'^ iterate of the solution x (and similarly the remainder r) are given by 
Xn = Xn-1 +'^«n^n-i- The global sums needed in the calculation of are done in double 
precision. Any residual errors, Sa^, can be absorbed into u and do not adversely affect 
the convergence rate as long as Un = ^(1 + dan) stays within the optimal range. Our tests 
also show that the calculation of can be done in single precision, and the convergence 
rate is not affected provided the remainder is refreshed somewhat more often depending 
on the quark mass. 

The only place that round-off errors due to use of 32 bit precision could affect the 
results is in the evaluation of the final convergence. For this purpose we make two checks. 
Firstly, we monitor the final value of |r^| and on each time-slice in addition to 

the global value which could be biased by time-slices closest to the source. We find that 
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the values at the time-shce farthest from the source are l^^^l/lx^l ~ 10~^ for the C quark 
and < 10"^"^ for other quarks. This means that, even for the heaviest quark, the effect 
of incomplete convergence is smaller than the statistical errors. Secondly, we double the 
number of inversion sweeps on randomly selected lattices. We find that the relative change 
in the value of hadronic 2-point correlation functions is < 10~^, which is negligible. 
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Fig. 1: Comparison of M^sit) for UiUi pion correlators with SL and WL sources. 



3. Fitting 

To illustrate some of the issues involved in fitting, we show, in Figs. representative 
results for the effective mass. We define this as Meff{t) = \n[C{t — l)/C{t)], plus corrections 
due to periodicity in the time direction. The source lies at t = in all plots. It is a general 
feature that the WL effective masses approach their asymptotic value from below, while 
those for SL and SS correlators approach from above. We consider in turn the problems 
we face and the solutions we have adopted. 
1. A major problem is that, at the 2a level, the convergence to the asymptotic value of 

Mgff is extremely slow. The SL and WL effective masses do come together eventually 
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Fig. 2: Comparison of Meff(t) for UiUi rho correlators with SL and WL sources. 



for pion and rho correlators, but for the nucleon the signal disappears into noise before 
convergence has occurred. Masses extracted from WL correlators are systematically 
lower than those from SL or SS correlators. For the pion and rho channels, this 
difference is ~ la, while for the nucleons it is ~ 2 — 3a. 

We think that this behavior is mainly due to limited statistics — a recent JLQCD study 
indicates that fluctuations at the 2a level are the rule rather than the exception. 
What we have done in practice is to average the masses with weightings {2*WL + SL + 
SS)/4. The equal weighting of wall and Wuppertal sources ensures that the resulting 
mass lies between the two sets of data on the effective mass plot. Errors are obtained 
using single-elimination jack-knife, performing a complete analysis (including forming 
the average masses) on each jack-knife sample. 
2. The previous problem is exacerbated by the fact that we have been unable to do stable 
fltting using the full correlation matrix. When we do so, some flts lead to clearly 
unreasonable values for masses, presumably because small errors in the correlation 
matrix are magnifled when this matrix is inverted to calculate x^- This problem is 
well known, and various remedies have been proposed [§] . For example, one can project 
against the eigenvectors of the correlation matrix having small eigenvalues, and then 
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Fig. 3: Comparison of Meff(t) for UiUiUi nucleon correlators with SL and WL sources. 



invert. Alternatively, one can reduce the range of times to which one fits. We have 
tried both schemes, but find that the resulting mass estimates are indistinguishable 
from those obtained fitting just with the diagonal elements of the correlation matrix 
( "uncorrelated fits"). Given that the schemes involve further subjective selections 
(e.g. how many eigenvectors to discard) we choose to use uncorrelated fits for our 
standard results. We do use the "correlated fits", however, to check that our fits are 
reasonable. Given that this is so, our jack-knife errors should be reliable. 

3. In a similar spirit, we have chosen to use single mass fits for our final answers. We 
have also made fits with two masses, but find that these fits require considerable 
tuning by hand, because our minimization routines tend to give the same result for 
the two masses for a few of the jack-knife samples. When we remove this problem, the 
resulting estimate for the lightest mass is consistent with that from the single mass 
fit, while the errors are 50% larger. Given the large number of correlators that we fit, 
it is impractical to do the tuning for each channel. 

4. Finally, we must choose a fit range. We choose the minimum time by inspecting the 
effective mass plots and deciding where the "plateau" region begins. The maximum 
time is then taken to be that at which the diagonal errors have roughly doubled 
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compared to the beginning of the plateau, or that at which the signal shows a clear 
break. For the pion correlator, where the errors do not grow with time, this means 
that we use all points beyond the minimum time. 

Our overriding criterion is to include as many time-slices as possible in the fits. We have 
not succeeded in developing a robust automated procedure that meets this objective when 
using two mass fits or incorporating the full covariance matrix. Since we analyze ~ 4000 
channels, it is not practical to tune the fitting of each by hand. Thus, for consistency, we 
use, for all channels, single mass fits keeping only the diagonal elements of the correlation 
matrix. We feel confident in both our central values and our error estimates, however, since 
in the many channels that we have fitted by hand, the results do not change significantly 
when we use the full correlation matrix and two mass fits. 



4. Quenching errors 

In the last few years it has been argued that, in the quenched approximation, the r( is 
a pseudo-Goldstone boson, and that r( loops give rise to chiral logarithms which make the 
quenched approximation singular in the chiral limit . The same methods ( "quenched 
chiral perturbation theory" ) have been used to estimate the errors due to quenching. It is 
important to test the predictions of this theory against numerical results. There is some 
supporting evidence — see Refs. ||Tl]] and [0 for recent reviews — but more work is needed. 



We collect here the results of quenched chiral perturbation theory relevant for this 
work. The expansion for the mass of a "pion" composed of quarks of mass m\ and m2 is 

Ml = c^(mi+m2)(l-51og[(mi + m2)]) + e^mJ + -- - , (4.1) 

with Ctt, e-K and b constants, and the ellipsis representing higher order terms in the chiral 
expansion. The term is shorthand both for analytic terms, i.e. those proportional to 
(mi-(-m2)^, rn\^ and m|, and for non-analytic chiral logarithms of general form log(mq). 
Such terms are present both in quenched and in full QCD, although the constants mul- 
tiplying them will be different in the two theories. In contrast, the b term is an artifact 
of quenching — it arises from 77' loops, and is divergent in the chiral limit. Evidence for 



this divergence has been found with staggered fermions [|T3|, but the effect is small, and 
becomes noticeable only at quark masses smaller than those we use. Thus when fitting the 
pion masses we ignore the 5 term. 

The predicted form for baryon masses is, schematically flj]. 



Mat = Oat -f b^^m}!'^ + h'^niq In(mg)] -|- c^niq 4- d]svr?J'^ ^r eis'rni H , (4.2) 



where b is the same constant as in Eq. ( [4.1J ), while — are additional constants. The 
expansion has the same form in full QCD, except that the b term, which again comes from 
rj' loops, is absent. As for the pion, we ignore the b term in almost all fits. It is small 
(because b is small), and, furthermore, is numerically indistinguishable from the higher 
order terms within our range of light quark masses. What we can test in some detail, 
however, is the expected form of the mg, m^^^, and terms. Here we benefit greatly 
from our use of non-degenerate quarks. 
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Similar comments apply to the vector meson masses. Although these have not been 
discussed in the quenched chiral perturbation theory literature, it is straightforward to 
extend the work done in QCD and arrive at the prediction 

Mp = ap + 5hpm]/'^ + Cp(mi + 1712) + dprr?J'^ + Cpm^ H . (4.3) 

Unlike for baryons, the detailed expressions for 6p, Cp and dp in terms of constants in the 
quenched chiral Lagrangian are not known. Nevertheless, we expect that the b term will 
be small, and again ignore it in our fitting. 
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Table 1. Meson masses from WL correlators at p = 0. 



5. Meson masses 

We give our results for meson masses in Tables |l|— We quote separately the results 
for the three source-sink combinations and the average (2 * WL + SL + SS)/4. There are 
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Table 2. Meson masses from SL correlators at p = 0. 



four columns for the pion in each table because we use Lorentz structures 75 = P and 
7475 = for both source and sink. Thus PP has a pseudoscalar interpolating field at 
both source and sink, while PA4 has an axial interpolating field at the sink. The four 
possibilities yield consistent results for masses, and we use the average of the four (inside 
our jack-knife loop) to give our best estimate of Mj^. 

States ai p = created by i/jiIj are labeled ao and those by i^'ji'y^il' are called ai. The 
signal in these channels is not good. We only present data for those mass combinations 
where there is a "plateau" over at least three time-slices. Even in the best case the signal 
dies out hj t = 14, so contamination from higher states is likely. We do not consider the 
data good enough to warrant further analysis. 

The tables show in detail how all the masses are systematically lower for the WL 
correlators than for the SL or SS correlators. The effect is, however, less than la for the 
pion and rho. As discussed above, we take the results of Table ^ as our best estimates. 

To extrapolate the hadron masses towards the chiral limit, and to test the forms 
predicted by quenched chiral perturbation theory, we need to choose a definition of quark 
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0.800(02) 
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0.792(03) 
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Table 3. Meson masses from SS correlators at p = 0. 



mass. We consider two possibilities. The first is the standard perturbative definition 

mf^ = Z^{1/2k-1/2kc), (5.1) 
while the second is non-perturbative [|16[ 



_ (0|A4(t)J(0)|0) Zp {mi,np + m2,np) , . 

" (0|P(t)J(0)|0) Za 2 ■ ^ ■ > 

Here P and A4 are local operators and J (which has the same Lorentz structure as either 
P or A4) is constructed from the smeared Wuppertal source propagators. The data for 
the two choices of J are consistent, so we use the average as the best estimate. We use 
tadpole improved renormalization constants, defined by []T7 



Z^ZIZa = Vl-3Ki/4/^ev/l-3K2/4/€c (l-0.316«^(g*)), (5-3) 



2 



'Z^ZlZp = - 3^i/4/^cVl - 3k2/4k, (1 + log W " 1-0335)) . 
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1.217(01) 


1.217(01) 


1.229(01) 


1.426(07) 


1.441(06) 


cs 


0.854(01) 


0.854(01) 


0.854(01) 


0.854(01) 


0.879(01) 




1.113(08) 


CUi 


0.814(01) 


0.814(01) 


0.814(01) 


0.814(01) 


0.842(01) 




1.086(10) 


CU2 


0.798(01) 


0.799(01) 


0.799(01) 


0.799(01) 


0.828(01) 




1.081(12) 


CUs 


0.789(01) 


0.791(02) 


0.790(01) 


0.790(01) 


0.819(02) 




1.084(15) 


ss 


0.422(01) 


0.422(01) 


0.422(01) 


0.422(01) 


0.505(01) 


0.727(12) 


0.763(08) 


SUi 


0.363(01) 


0.363(01) 


0.364(01) 


0.363(01) 


0.464(01) 


0.700(16) 


0.730(08) 


SU2 


0.339(01) 


0.339(01) 


0.339(01) 


0.338(01) 


0.448(02) 


0.718(21) 


0.720(11) 


SUs 


0.323(01) 


0.323(01) 


0.323(01) 


0.323(01) 


0.439(02) 




0.720(12) 


UlUi 


0.296(01) 


0.296(01) 


0.297(01) 


0.296(01) 


0.422(02) 


0.699(31) 


0.694(10) 


UIU2 


0.267(01) 


0.266(01) 


0.267(01) 


0.266(01) 


0.405(02) 




0.683(12) 


UIU3 


0.247(01) 


0.247(01) 


0.248(01) 


0.246(01) 


0.394(03) 




0.682(15) 


U2U2 


0.234(01) 


0.233(01) 


0.234(01) 


0.233(01) 


0.387(03) 




0.672(14) 


U2US 


0.211(01) 


0.210(01) 


0.211(01) 


0.210(01) 


0.373(03) 




0.669(17) 


UsUs 


0.185(01) 


0.184(01) 


0.186(01) 


0.184(01) 


0.361(05) 




0.664(20) 



Table 4. Best estimates of meson masses (given by the average (2 * WL + SL + SS)/4 as 
explained in the text) at p — 0. 



1 2 

Z^' are wavefunction renormalizations for the two quarks coupling to the bilinears, which 
have hopping parameters ki^2 respectively. The Z^'s cancel in the ratio Zp/Za- p is the 
scale at which we match to the continuum renormalization scheme (here MS), while q* 
is a typical momentum in the 1-loop integral. We choose /i = 1/a, and assume q* = 1/a. 
The masses can be run to other scales using the two-loop anomalous dimension relation 

m 

1 + T7rr2 ( 77^2 ) ■ (5-4) 



m(^) \gH^^)J V 167r2 ^ 2/3^ 

We mostly quote results for the scale 1/a, but also give some results for the scale 2 GeV, 
to allow comparison with other work. 

It is important to note that, although we do not label m„p with MS explicitly, both 
definitions yield estimates for the continuum mass in the MS scheme, and should agree in 
the continuum limit. We give our results for both definitions in Table |^, and find a sub- 
stantial difference between them. Most likely this is due to the failure of the perturbative 
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expression for Zp [IT]. Details of this analysis will be presented elsewhere. 

This difference is not, however, important for our study, since the two definitions 
of mass turn out to proportional to very high accuracy for the 6 cases {UilJj\. This is 
shown in Fig. ^, in which we plot the average m^p of the quark and anti-quark against 
the average value of 1/2k. The fit is given by rrin-p = 0.867(4) (1/2k) — 2.76(1). Thus, 
when extrapolating, it does not matter which definition of mass one uses. We choose to 
use vrin-p for most purposes in this paper — and we will drop the subscript henceforth. The 
perturbative mass, when it appears, will be distinguished by the superscript MS. 

We now consider extrapolations to the chiral limit. Figures ^ and show, respectively. 
Ml and Mp function of the average m^p- Only data from the combinations {UilJj\ 
(the six points with the smallest masses) and {SS, SUj} (the rightmost four points) are 
included. 

We first consider linear fits to the data. These work well for the {UiUj} points, and 
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are the fits shown in the figures. The fits give 

= + c^(mi +7712) = 0.0013(5) + 2. 296(ll)(mi +m2)/2, 

(5.5) 

Mp = ap + Cp{mi+m2) = 0.327(6) + 2.54(14)(mi + m2)/2 . 

One interesting feature is that the pion intercept is slightly inconsistent with zero, whereas 
in the continuum limit of full QCD the intercept should vanish. This discrepancy could 
be due to discretization errors: for to be proportional to rrinp it must be true that 
(OIA4IP) oc (see Eq. (|5.2|) ) which is guaranteed only in the continuum limit. A violation 
of this proportionality at 0{a) would lead to an additional term in Eq. (|5.5| ) proportional 
to yjarrinp, where a is the lattice spacing. This could give rise to a non-zero intercept. 
Other possible culprits are quenched chiral logarithms — the terms proportional to 5 in 
Eq. (O) — and finite size effects. We expect the latter to be insignificant, however, since 



Mt^L > 6. Our data is not good enough to investigate this effect further. As a result of 
the non-zero intercept, the chiral limit (defined by = 0) is at m^J**'^"^ = -0.00058(21). 
In the following, when quoting quark masses in physical units we use rrinp — m^p*'^"^, i.e. 
we offset the zero of the mass scale. This offset is only significant for m. When we use 
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masses in lattice units, as we do in all our fits and extrapolations, we do not include this 
offset. 

Using our linear fits we determine frinp by extrapolating to the point where M^/M^ = 
(137/768)^, finding Winp = 0.00093(21). The resulting value of MpU gives us an estimate 
of the scale 

a-i (Mp) = 2.330(41) GeV . (5.6) 

Using this value, and including the offset in the mass scale, we find m(l/a) — 3.51(6) MeV 
(m(2 GeV) = 3.60(5)). To facilitate comparison with earlier calculations it is useful to 
give results for the critical and light quark hopping parameters. Using linear fits versus 
1/k we find 

= 0.157131(9), = 0.157046(9) . (5.7) 

We note that linear fits work well for the subsample {SS, SUj} (and for {CUj}, not shown 
in the figures) as well, and we use them below when calculating the masses of strange and 
charm mesons. 

The figures show that linear fits are inadequate for the combined {UiUj} and 
{SS, SUj} data. This is particularly striking for Mp, where there is a definite nega- 
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tive curvature, which, if ignored, could lead to an overestimate of ap. Thus we have tried 
fits to the combined data involving higher order terms in the chiral expansion. A quadratic 
fit to works well, but has little impact on extrapolations, increasing a^r by about la. 
For the p, motivated by chiral perturbation theory, we have used two forms of higher order 
terms 

Mp = ap + Cp{mi + 7712) + dp{mi + m^)'^^'^ , 

(5.8) 

-^p = ftp + Cp(mi + 7712) + ep(mi + m2)^ . 

These fits, when extrapolated to m (calculated self-consistently), yield smaller values of 
Mp(m) and thus larger scales than the linear fit: 

a-^(M„) = 2.365(48) GeV (m^/^ fit) , 

(5.9) 

a-i(Mp) = 2.344(42) GeV {n? fit) , 



The difference from the result of the linear fit is, however, within the statistical errors. Our 
data is good enough to see the curvature, but not precise enough to study it in detail. In 
particular, we cannot distinguish between the two forms in Eq. ( p.8| ), both doing a good 
job of fitting the combined {UiUj\ and {5* 5", S\Jj\ data sets. In view of this, we chose to 
use the scale from the linear fit for most of our subsequent analysis. For comparison, we 
note that the scale we find from /tt is 2265(57) MeV [T^, while the NRQCD collaboration 



reports a value 2.4(1) GeV from a mean of charmonium and T spectroscopy [^. Thus, 
different estimates based on meson correlators are consistent. In the remainder of the paper 
we shall take the scale from Mp and assign a 3% systematic uncertainty to cover the spread 
in obtained from different mesonic observables and different types of extrapolations. 

Using the linear fits defined in Eq. ( p.5|) we determine the lattice strange quark mass 
by first extrapolating M^/M^, Mk* /Mp, or M^jMp to m, and then linearly interpolating 
between \J\ and 5" until these quantities match their physical values. (We use Mk = 
495 MeV, Mk* = 894 MeV, and M^ = 1019 MeV.) The results are given in Table | in 
terms of n and the two definitions of quark mass discussed above, except that we have here 
run the masses to a scale of 2 GeV to facilitate comparison with other work. 



Qty. 


Kg 


ms,np{l/a) 


ms,np{'^ GeV) 


1 

mf^(2 GeV) 


Ml/Ml 


0.15503(7) 


87(2) 


89(2) 


129(2) 


Mk*/Mp 


0.15479(19) 


98(7) 


100(7) 


145(9) 


M^/Mp 


0.15464(17) 


104(6) 


106(6) 


154(8) 



Table 5. Estimates of the strange quark mass obtained by matching different quantities 
with their physical values. Results are given for (1) Kg (2) the lattice non- 
perturbative mass mnp/a evaluated at 1/a and 2 GeV, and (3) the lattice per- 
turbative mass mf^^ evaluated at 2 GeV. The latter three masses are in MeV. 
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The three ratios lead to significantly different results for nis, presumably because of a 
combination of quenching and discretization errors. Using M'^/M^ to fix nis implies = 
25m, since we use the lowest order chiral expansion to fit the data. On the other hand, 
the estimates using Mfj^/Mp and MK*/Mp are not constrained by the chiral expansion, 
and give nis/rn ~ 30, in surprisingly good agreement with the next-to-leading chiral result 
pi}] . In this paper we quote all results using ms(M^). 

Recently Lacock and Michael [|2| suggested using the dimensionless quantity Jy = 
MydMv /dM^ to test the quenched approximation. Using a linear fit we find 

Jk* = 0.41(1), 

5.10 

Jp = 0.36(1), 

to be compared to the experimental values of ~ 0.48 and ~ 0.41 respectively. Again, the 
discrepancy is presumably due to a combination of discretization and quenching errors. 



6. 0{ma) discretization errors 

At /3 = 6, the charm mass is such that nica ~ 1, so there are potentially large 0{mca) 
discretization errors in all quantities involving charm quarks. These are in addition to 
errors of ©(Aqcd^): which are common to all quantities. One effect of 0(a) errors is that 
there is a difference between the "static" mass Mi = E{'p = 0), and the "kinetic" mass 
M2 = {d"^ E / dp^\p=Q)~^ . Here the energy is determined from the rate of exponential decay 
of the correlator, C{t) oc exp(— i?(p)i). Mi and M2 agree in the continuum limit, whatever 
the mass of the state. 

To evaluate M2, we have tested four forms of dispersion relation 

{A) E^ =p'^ + M^ ^ M2 = M ; 

(S) sinh^ E = sin^ p + sinh^ M ^ M2 = - sinh 2M ; 

E v M (6-1) 

(C) sinh^ — = sin^ - + sinh^ — M2 = sinh M ; 

2 2 2 

(D) smh^ — = sm^ - + smh^ — M2 = 2 smh — . 

^ ' 4: 4: 4 2 

(A) is the continuum relation, while (B), (C) and (D) are lattice forms following from 
different choices of lattice action. In particular, (C) follows from the nearest neighbor 
symmetric difference discretization of the action for a scalar. Our results for E{p) are 
collected in Table ^, and in Fig. ^ we show how the various dispersion relations fare for 
the CUi meson. For all heavy-heavy and heavy-light mesons, our results turn out to be 
consistent with dispersion relation (C), but not with the other forms. 
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p2 = i 


p2 = 2 


p2 = 3 


= 4 


p2 = l 


p2 = 2 


p2 = 3 


= 4 


cc 


1.231(01) 


1.244(01) 


1.255(02) 


1.268(01) 


1.243(01) 


1.255(01) 


1.267(02) 


1.278(02) 


cs 


0.874(01) 


0.893(02) 


0.913(03) 


0.930(02) 


0.899(01) 


0.917(02) 


0.937(03) 


0.951(03) 


CUi 


0.836(02) 


0.855(02) 


0.876(03) 


0.895(02) 


0.863(02) 


0.881(02) 


0.901(03) 


0.917(04) 


CU2 


0.821(02) 


0.841(02) 


0.861(04) 


0.881(03) 


0.849(02) 


0.868(03) 


0.888(04) 


0.905(04) 


CUs 


0.814(03) 


0.833(03) 


0.853(04) 


0.873(03) 


0.841(03) 


0.860(04) 


0.882(04) 


0.898(05) 


ss 


0.466(01) 


0.504(02) 


0.547(04) 


0.578(04) 


0.541(02) 


0.575(03) 


0.614(05) 


0.642(04) 


SUi 


0.414(02) 


0.456(03) 


0.505(06) 


0.539(04) 


0.503(03) 


0.540(03) 


0.583(06) 


0.613(05) 


SU2 


0.393(02) 


0.438(03) 


0.487(07) 


0.522(05) 


0.489(03) 


0.528(04) 


0.571(07) 


0.602(06) 


SUs 


0.380(02) 


0.427(03) 


0.474(08) 


0.511(06) 


0.482(04) 


0.522(05) 


0.565(09) 


0.596(07) 




0.357(02) 


0.406(03) 


0.460(08) 


0.497(06) 


0.465(04) 


0.508(05) 


0.553(08) 


0.583(06) 


UIU2 


0.333(02) 


0.386(03) 


0.438(09) 


0.480(07) 


0.452(05) 


0.496(06) 


0.541(10) 


0.572(07) 




0.317(02) 


0.372(04) 


0.423(11) 


0.469(08) 


0.445(06) 


0.490(07) 


0.534(12) 


0.566(08) 


U2US 


0.307(02) 


0.363(04) 


0.416(11) 


0.464(09) 


0.439(06) 


0.484(08) 


0.529(12) 


0.561(09) 


U2US 


0.290(03) 


0.349(05) 


0.399(13) 


0.454(11) 


0.432(08) 


0.478(10) 


0.520(15) 


0.554(10) 




0.272(03) 


0.333(05) 


0.379(16) 


0.445(13) 


0.424(10) 


0.471(13) 


0.509(19) 


0.546(13) 



Table 6. Pseudoscalar (first set) and vector (second set) meson energies as a function of 
momentum. The data are tlie average of SL and SS estimates. 



Using this result, we have that M2 = sinhMi. In Table we compare Mi and M2 for 
the pseudoscalar and vector mesons. The difference is tiny for the smallest quark masses, 
but substantial for charmed mesons. Our results for charmed meson masses (obtained by 
linear extrapolation in the light quark mass and with mc chosen to be k = 0.135) are 
given in Table |^. There is a significant difference between the estimates using Mi and 
M2: for the D mesons the difference is ~ 10%, while for the charmonium system it is 
25 — 30%. This suggests that 0{ma) corrections in the Wilson action are already large 
in the charm region. This is also apparent from the mass-splittings — the spin-spin and 
spin-orbit interactions are underestimated by the Wilson action, as has been previously 
observed [ p3| . 

7. Baryon masses 

Our baryon mass analysis is based on two overlapping data sets. For the complete 
data set (170 configurations) we have results for only WL and SL correlators, and for only 
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Fig. 7: Test of four lattice dispersion relations (see Eq. 1\) ) using the 
lattice pion data at various momenta for the case CUi . The data 
favor the nearest-neighbor symmetric-difference relativistic disper- 
sion relation sinh^(i?/2) — sin^(p/2) = sinh^(M/2), as shown by the 
square symbols. 



a subset of possible quark combinations. On the last 110 configurations we also calculate 
SS correlators, and use all degenerate and non-degenerate combinations made up of S and 
Ui quarks. Our results for the masses are given in Tables ^, |TD|, and ^ The sample size is 
indicated by the subscript in the table headings. In the following, some analyses are only 
possible on the smaller sample, and we label such results by an asterisk next to the error 
estimate. 



7.1. Correlators 

It is necessary to introduce some notation to explain the correlation functions we 
calculate. We adapt that used in quenched chiral perturbation theory for baryons jl^ . 



Assume that we have three flavors of quark, labeled by w, d and s. To create spin-1/2 
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State 


Ml 




Ml 


Mo 






CC 


L.ZL 1 [\J) 




i.zzy(^i J 


1 c;fi/i/'9'i 

i.0D4(^Z J 


0.4569(2) 


0.6563(2) 


cs 


U. 504(^1 J 




u.o / y(^i J 


u.yyy (^z j 


0.2625(2) 


0.3822(2) 


CUi 


U. 514(^1 ) 




U.o4Z(^i J 


u.y4 < (^zj 


0.2419(3) 


0.3557(2) 


CU2 


u. I yy^^i ) 




u.ozo^^i J 


u.yz / (^o ) 


0.2339(3) 


0.3453(2) 




U. ItfUyl ) 


U.o ( ( \Oj 


u.oiy (^z J 


U.y iD(^4 J 


0.2291(4) 


0.3388(2) 


ss 


n Aoo(^ ^ 




U.OUO(^l J 


n f^98('9'\ 
U.OZo(^Z J 


0.0756(2) 


0.1081(2) 


SUi 


n qrqM ^ 
U.oDo(^l ) 


n Q70M ^ 
u.o ( 


U. 404(^1 J 


n /I81 ('Q'l 


0.0565(2) 


0.0816(2) 


SU2 




A Q /I ft /I \ 


A /I /I Q /'0'\ 

U.44o(^Z J 


U.4D0(^o j 


0.0490(2) 


0.0711(2) 






U.ooUl^i J 


u.4oy(^z J 


U. 404(^4 J 


0.0445(2) 


0.0647(2) 


U^U^ 








U.40U1^0 J 


0.0377(2) 


0.0550(2) 


UIU2 


0.267(1) 


0.271(1) 


0.405(2) 


0.417(4) 


0.0304(2) 


0.0446(2) 


UlUs 


0.247(1) 


0.250(1) 


0.394(3) 


0.406(5) 


0.0259(2) 


0.0382(2) 


U2U2 


0.234(1) 


0.236(1) 


0.387(3) 


0.398(5) 


0.0232(2) 


0.0342(2) 


U2US 


0.211(1) 


0.213(1) 


0.373(3) 


0.385(7) 


0.0187(2) 


0.0277(2) 


UsUs 


0.185(1) 


0.186(1) 


0.361(5) 


0.370(9) 


0.0143(2) 


0.0213(2) 



Table 7. Comparison of Mi = E{j) = 0) and M2 = sinhMi. We also give the values of 
the average non-perturbative rrinp and perturbative definitions of the quark mass 
in the MS scheme as defined in Eqs. ( |5.2| ) and ( \5.3i ). Both are in lattice units, 
evaluated at the scale 1/a. 



baryons we can use the interpolating operators 

o^,,)k = ii^l.c^^M^cM""' , (7.1) 

where a, b and c label color, while i, j and k label flavor. It is simple to show that 
^{ij)k = ~0{ji)k^ SO that there are only nine independent operators — eight SU{?>) octets 
and the singlet Xlij/c ^^''^^ijk- One way to project against the singlet is to form 

Bijk = 0(^ij)k + C>{ik)j = Bikj ■ (7.2) 
There are eight independent Sijfc's, the relation to the usual states being exemplifled by 

(7.3) 

id ■ 



The overall factor in these equations is arbitrary, while the relative normalization is flxed 
by SU{3) symmetry. 



V2p = 


Bduu ^Bfiud 




BsuU ^BfiUS 


eo = 


Buds + Bdsu = —B 


y3A° = 


Buds Bdsu • 
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D meson masses in MeV 




IVl \ 


IVI2 


Hi 'V 'KIT 

HiJbpi. 






iyyU(^o4 J 


ioby 


ivij:)* 




zUobi^oo J 


zUUo 




1 QCid. ( Q^^ 

ioybl^oU ) 


01 10/' Qo^ 
ziiz(^oz j 


iyby 




1 ni A ( o^^^ 


Zii) ( (z (^J 


1969 


MD*{ms[MK)) 


1 c\(i 1 / o 1 ^ 


zz0i(o4J 


ziiO : 


A/T I I A/T W 


iy ( o(^z ( ) 


zzz4(^zyj 


01 1 n9 

ziiU ; 


A/T ('1 c ^ 




3590(64 J 


zyoO 


A/T ('3 C ^ 


ZOD0(0i ) 


364o(65 J 


309 ( 


A/T f'i r> ^ 


ooZ4(^DU J 


40 / Z(^oD J 


O4i0 




3357(60) 


4646(86) 


3510 


AM{^Si - ^So) 


29(1) 


53(2) 


117 


AM(3Pi - 3Po) 


33(9) 


74(18) 


95 


AM(3Po - ^^i) 


459(17) 


929(36) 


318 



Table 8. A comparison of lattice estimates of D meson and charmonium masses with the 
experimental data. We show results for Mi and M2 and for the two different 
ways of setting described in the text. 



All spin-1/2 baryon correlators are built out of the two contractions shown in Fig. ^ 
The notation {DU)S = {UD)S corresponds to quarks of flavors U and D contracted into 
a closed loop, while the propagator for 5" carries the spin quantum numbers of the baryon. 
The notation [DUS) corresponds to a single ordered contraction of the three quarks. We 
consider two types of correlator, "E-like" and "A-like" . The former is exemplified by that 

of the E° 

S{UD} = S{DU} = {B,^d{x) BsudiO)) 

(7.4) 

= {US)D + {DS)U + (USD) + (DSU) . 

(This equation defines our sign conventions for the contractions.) The proton, neutron, 
E+, E~, and E~ correlators are also of this type: they are, respectively, D{UU}, 
U{DD}, S{UU}, S{DD}, U{SS} and D{SS}. The second type of correlator is that of 
the AO 

spn, ^ sm, . ( V3) ms^n . iosw . Mums - iuso, - (z,.., 

+2{SUD) + 2{SDU) + 2{UDS) + 2{DUS)] . 
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When rriu rud, there is also a non- vanishing A° — E*^ cross correlator, but we have not 
found this to give useful results. 

We have calculated the two types of spin- 1/2 baryon correlator for all independent 
mass combinations involving Ui and 5" quarks. The results are given in Tables ^ and p!0| 
respectively. Masses from correlators of the form A [AS] and A{AB} are also included 
even though they are not independent — the correlators are related by 

A[AB] = - {3B[AA] + B{AA}) , 

f (7.6) 
A{AB} = - {B[AA] + 3B{AA}) . 

One can think of the results for the A{BB} and masses as being those for the 

S and A, respectively, with nis = rriA and = = m^. Unlike in the real world, 
there is nothing to stop these two masses being the same, i.e. = ms, in the quenched 
approximation. Note, however, that in this case the S and A are also degenerate, i.e. 
M{A{AA}) = M{A[AA]). Indeed, the contractions in the two cases are identical. 

The interpretation of the results for the completely non-degenerate correlators, A[i?C] 
and A{BC}, is more complicated. Because isospin is broken, the E*^- and A-like states 
mix, with both correlators containing contributions from both physical states. Let M_|_ 
and M_ be the masses of the heavier and lighter states, respectively, and dM the mass 
difference. At long times, the effective mass for both correlators will asymptote to M_. 
However, at times short compared to the inverse mass difference, i.e. dMt^a.x ^ 1, there 
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5'-^iio 


5''S'iio 


AViio 




0.786(03) 


0.791(04) 


0.793(04) 


0.789(02) 


S[SUi] 


0.734(04) 


0.740(04) 


0.742(04) 


0.738(02) 


S[SU2] 


0.712(04) 


0.720(04) 


0.721(04) 


0.716(03) 


S[SU3] 


0.698(04) 


0.707(05) 


0.708(04) 


0.703(03) 


S[UiUi] 


0.678(01) 


0.687(01) 


0.688(01) 


0.683(03) 


S[UiU2] 


0.655(04) 


0.665(05) 


0.665(04) 


0.660(03) 


S[UiUs] 


0.640(05) 


0.652(05) 


0.651(05) 


0.646(03) 


S[U2U2] 


0.630(05) 


0.642(05) 


0.642(05) 


0.636(03) 


S[U2Us] 


0.614(06) 


0.627(06) 


0.627(05) 


0.620(04) 


S[UsUs] 


0.597(07) 


0.611(07) 


0.610(06) 


0.604(04) 


Ui[SS] 


0.740(04) 


0.747(04) 


0.749(04) 


0.744(03) 


Ui[SUi] 


0.687(05) 


0.696(05) 


0.697(04) 


0.691(03) 


Ui[SU2] 


0.663(04) 


0.675(05) 


0.675(05) 


0.669(03) 


Ui[SUs] 


0.648(05) 


0.661(05) 


0.661(05) 


0.655(03) 


Ui[UiUi] 


0.630(05) 


0.643(05) 


0.642(05) 


0.636(03) 


Ui[UiU2] 


0.605(05) 


0.620(06) 


0.619(05) 


0.612(03) 


UliUiUs] 


0.589(06) 


0.605(06) 


0.604(06) 


0.597(04) 


Ul[U2U2] 


0.579(06) 


0.596(06) 


0.595(06) 


0.587(04) 


Ui[U2U3] 


0.562(07) 


0.580(07) 


0.579(06) 


0.571(04) 


Ui[UsUs] 


0.544(08) 


0.563(09) 


0.562(07) 


0.553(05) 


U2[SS] 


0.723(04) 


0.731(05) 


0.732(04) 


0.727(03) 


U2[SU-l] 


0.668(05) 


0.679(05) 


0.679(05) 


0.671(03) 


U2[SU2] 


0.643(05) 


0.657(06) 


0.657(05) 


0.650(03) 


U2[SUs] 


0.627(06) 


0.643(06) 


0.643(06) 


0.635(04) 


U2[UiUi] 


0.610(06) 


0.626(06) 


0.625(05) 


0.617(04) 


U2[UlU2] 


0.584(06) 


0.602(07) 


0.601(06) 


0.593(04) 


U2[UiU^] 


0.567(07) 


0.587(08) 


0.586(07) 


0.577(05) 


U2[U2U2] 


0.557(07) 


0.578(08) 


0.577(07) 


0.567(05) 


U2[U2U3] 


0.539(08) 


0.561(09) 


0.560(08) 


0.550(05) 


U2[UsUs] 


0.520(10) 


0.543(11) 


0.543(09) 


0.532(06) 


Us[SS] 


0.712(05) 


0.722(05) 


0.723(05) 


0.717(03) 


Us[SUi] 


0.656(06) 


0.670(06) 


0.669(05) 


0.663(04) 


U3[SU2] 


0.631(06) 


0.647(07) 


0.646(06) 


0.639(04) 


Us[SUs] 


0.614(07) 


0.632(08) 


0.631(07) 


0.623(04) 


Us[U,Ui] 


0.598(06) 


0.616(07) 


0.615(06) 


0.607(04) 


Us[UiU2] 


0.572(07) 


0.592(08) 


0.590(07) 


0.581(05) 


U3[UiUs] 


0.554(08) 


0.576(09) 


0.574(08) 


0.564(05) 


U3[U2U2] 


0.544(08) 


0.567(09) 


0.565(08) 


0.555(05) 




0.525(09) 


0.550(11) 


0.519(09) 


0.537(06) 




0.506(11) 


0.531(13) 


0.531(10) 


0.518(07) 



Table 9. Mass estimates for A-like baryons. 
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will be an approximate plateau at a value which is a weighted average of the two masses. 
To see this, we pick the A correlator and write it as 



CA{t) = Ae-'^'-*(cos'^ + sin'^e-'^'^'*) + ■ ■ ■ , (7.7) 

where tan^ is the ratio of the amplitudes to create the two mixed states, and the ellipsis 
represents excited states. The effective mass is 

m(A)eff(t) = _ ^^'^^^(^) ^ M- + siii^e 6M{1 + 0{6Mt)) ^ cos'^O M- + sin^e M+ . (7.8) 

Thus the effective mass is almost constant, and given our errors, we cannot distinguish it 
from a plateau. We discuss below the interpretation of the resulting "mass" . 

We have also calculated the masses of spin-3/2 baryons. Here there is only one type of 
operator, which is completely symmetric in flavor. Our results for all 20 mass combinations 



built out of Ui and 5" quarks are given in Table |TI. 

We analyze the baryon masses using the chiral expansion. We first consider the spin- 
1/2 baryons, returning to spin-3/2 baryons later in this section. We couch the discussion 
in terms similar to those used in full QCD. If we keep only constants and terms linear in 

1/2 

quark masses, (and drop non-analytic terms of the form 5mq as discussed in section H), 



then it is straightforward to show, using quenched chiral perturbation theory [0, that 



M(E+) = M{S{UU}) = Mo + 4Fm^ + 2{F - D] 



m 



S 5 



2D D (7-9) 

M(A) = M{S[UU]) = Mo + 4{F - — )m„ + 2{F + -)m, , 

Here Mq is spin-1/2 baryon mass in the chiral limit, and F and D are the usual reduced 
matrix elements of scalar densities. Note that there is no dependence on m^, since the 
d quark does not enter the correlators. These results can also be obtained using mass 
perturbation theory in full QCD, and then deleting terms proportional to ruu + ma + rus 
which correspond to contributions from internal quark loops. We stress that these formulae 
apply to all the states we consider — for example, the proton mass is M{D{UU}), and is 
obtained by replacing with in the formula for M(E"'"). 

At this order in the chiral expansion, it is simple to extend the results to baryons com- 
posed of three non-degenerate quarks. The mass matrix in the (E°, A) — [S{UD}, S\UD]) 
basis is 



a ^\ _ f Mo + 4Fm + 2{F - D)ms ^{niu - ma) 

1 ^(mu-md) Mo + 4{F-^)m + 2{F + ^)m, 



(7.10) 



Diagonalizing this matrix gives the eigenvalues M± with a mixing angle 9. If we assume 
that the same mixing angle applies for the interpolating fields, which we think is true up 
to corrections of 0(mq^^), then 9 is the angle appearing in our previous expressions for 
the A correlator (Eq. ( [7.7| )). The "short-time effective mass", Eq. ( |7.8D , is then 

M(A)eff ^ cos^^M. +sin2^M+ = /3. (7.11) 
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WLiio 


SLiio 


SSiio 


AViio 


WLi7o 


SLno 


AVno 


S{SS} 


0.786(03) 


0.791(04) 


0.793(04) 


0.789(02) 


0.786(03) 


0.794(03) 


0.790(02) 


S{SUi} 


0.738(04) 


0.745(04) 


0.746(04) 


0.742(03) 








S{SU2} 


0.718(04) 


0.727(04) 


0.728(04) 


0.723(03) 








S{SU3} 


0.706(04) 


0.716(05) 


0.717(05) 


0.711(03) 








S{UiUi} 


0.689(04) 


0.699(05) 


0.700(04) 


0.695(03) 


0.689(04) 


0.702(04) 


0.695(03) 


S{UiU2} 


0.669(05) 


0.682(05) 


0.682(05) 


0.676(03) 








S{UiUs} 


0.657(05) 


0.672(06) 


0.671(05) 


0.664(04) 








S{U2U2} 


0.649(06) 


0.665(06) 


0.664(06) 


0.657(04) 


0.649(04) 


0.667(05) 


0.658(03) 


S{U2U3} 


0.636(07) 


0.654(07) 


0.654(06) 


0.645(04) 








S{UsUs} 


0.622(08) 


0.644(09) 


0.643(08) 


0.633(05) 


0.625(06) 


0.647(07) 


0.636(04) 


Ui{SS} 


0.732(04) 


0.738(04) 


0.739(04) 


0.736(02) 


0.732(03) 


0.744(04) 


0.738(02) 


Ui{SUi} 


0.681(05) 


0.690(04) 


0.690(04) 


0.686(03) 








Ui{SU2} 


0.659(04) 


0.671(05) 


0.671(04) 


0.665(03) 








UiiSU^} 


0.646(05) 


0.659(05) 


0.658(05) 


0.652(03) 








Ui{lhUi} 


0.630(05) 


0.613(05) 


0.612(05) 


0.636(03) 


0.630(05) 


0.615(05) 


0.638(03) 


Ui{UiU2} 


0.608(05) 


0.624(06) 


0.623(05) 


0.616(04) 










0.594(06) 


0.612(07) 


0.611(06) 


0.603(04) 








Ul{U2U2} 


0.586(06) 


0.605(07) 


0.604(06) 


0.595(04) 


0.588(05) 


0.607(06) 


0.598(04) 


Ui{U2Us} 


0.572(07) 


0.593(08) 


0.592(07) 


0.582(05) 








Ui{UsUs} 


0.557(09) 


0.581(10) 


0.580(08) 


0.569(06) 


0.562(07) 


0.583(08) 


0.573(05) 


U2{SS} 


0.710(04) 


0.717(04) 


0.718(04) 


0.714(03) 


0.708(04) 


0.720(04) 


0.714(03) 


U2{SUi} 


0.657(05) 


0.667(05) 


0.667(04) 


0.662(03) 








U2{SU2} 


0.634(05) 


0.646(05) 


0.646(05) 


0.640(03) 








U2{SUs} 


0.619(06) 


0.633(06) 


0.633(05) 


0.626(04) 








U2{UiUi} 


0.603(05) 


0.618(06) 


0.618(05) 


0.611(03) 


0.605(06) 


0.621(06) 


0.613(03) 


U2{UiU2} 


0.580(06) 


0.598(07) 


0.597(06) 


0.589(04) 








U2{UiUs} 


0.565(07) 


0.585(08) 


0.584(06) 


0.575(04) 








U2{U2U2} 


0.557(07) 


0.578(08) 


0.577(07) 


0.567(05) 


0.561(06) 


0.581(07) 


0.571(04) 


U2{U2Us} 


0.542(08) 


0.565(09) 


0.564(08) 


0.553(05) 








U2{U3Us} 


0.526(10) 


0.552(11) 


0.551(09) 


0.538(06) 


0.534(08) 


0.555(09) 


0.544(06) 


UsiSS} 


0.696(05) 


0.704(05) 


0.704(04) 


0.700(03) 


0.694(04) 


0.708(05) 


0.701(03) 


U3{SUi} 


0.612(06) 


0.653(05) 


0.652(05) 


0.617(03) 








Us{SU2} 


0.617(06) 


0.631(06) 


0.630(05) 


0.624(04) 








Us{SUs} 


0.602(07) 


0.616(07) 


0.616(06) 


0.609(04) 








Us{UiUi} 


0.586(06) 


0.602(06) 


0.602(06) 


0.594(04) 


0.589(06) 


0.606(06) 


0.597(04) 


U3{UiU2} 


0.563(07) 


0.581(08) 


0.580(06) 


0.572(04) 










0.547(08) 


0.567(09) 


0.566(07) 


0.557(05) 








U3{U2U2} 


0.539(08) 


0.560(09) 


0.559(07) 


0.549(05) 


0.543(08) 


0.564(08) 


0.553(05) 


U3{U2U3} 


0.523(09) 


0.545(11) 


0.545(09) 


0.534(06) 










0.506(11) 


0.531(13) 


0.531(10) 


0.518(07) 


0.515(10) 


0.536(10) 


0.525(07) 



Table 10. Mass estimates for T,-like baryons. 
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SLiio 


SSiio 




WLiro 


SLiYo 


AVuo 


{SSS} 


0.832(05) 


0.843(06) 


0.843(07) 


0.837(04) 


0.831(04) 


0.845(05) 


0.838(03) 


{SSUi} 


0.788(06) 


0.803(07) 


0.802(07) 


0.795(05) 








{SSU2} 


0.770(07) 


0.789(07) 


0.786(08) 


0.779(06) 








{SSUs} 


0.760(08) 


0.781(09) 


0.778(08) 


0.770(06) 








{SUiUi} 


0.744(08) 


0.764(08) 


0.761(08) 


0.753(06) 








{SUiU2} 


0.725(10) 


0.750(09) 


0.747(08) 


0.736(07) 










0.710(11) 


0.741(11) 


0.738(09) 


0.725(08) 








{SU2U2} 


0.705(10) 


0.736(11) 


0.732(09) 


0.719(07) 








{SU2US} 


0.690(11) 


0.727(13) 


0.724(10) 


0.708(08) 








{SUsUs} 


0.674(09) 


0.718(15) 


0.716(11) 


0.696(08) 








{UlUiUi} 


0.697(08) 


0.726(10) 


0.721(09) 


0.710(06) 


0.699(08) 


0.721(09) 


0.710(06) 


{UIUIU2} 


0.677(09) 


0.712(12) 


0.706(10) 


0.693(07) 








{UlUiUs} 


0.663(11) 


0.704(14) 


0.699(11) 


0.682(08) 








{UIU2U2} 


0.656(09) 


0.699(14) 


0.693(11) 


0.676(08) 








{U1U2US} 


0.642(10) 


0.690(17) 


0.685(12) 


0.665(09) 








{UlUsUs} 


0.629(09) 


0.682(21) 


0.678(13) 


0.654(09) 








{U2U2U2} 


0.635(08) 


0.686(18) 


0.680(12) 


0.659(08) 


0.644(11) 


0.683(12) 


0.664(08) 


{U2U2US} 


0.622(10) 


0.678(22) 


0.673(14) 


0.649(10) 








{U2UsUs} 


0.609(11) 


0.670(27) 


0.667(16) 


0.639(11) 








{UsUsUs} 


0.596(13) 


0.661(33) 


0.662(18) 


0.628(14) 


0.612(13) 


0.660(22) 


0.636(13) 



Table 11. Mass estimates for the decuplet baryons. 



A similar argument shows that 

M(S)eff ^ sin^ OM- + cos^ eM+ = a . (7.12) 

Thus we find the surprising result that the short-time effective masses are insensitive to 
the isospin breaking term 7. Furthermore, the expressions for a and /3 are exactly the 
same as the formulae applicable when isospin is unbroken, Eqs. (^]D, except that is 
replaced by the average mass, m. 
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We have not extended this analysis to higher order in the chiral expansion. Neverthe- 
less, we use it as motivation for including our results for baryons composed of completely 
non-degenerate quarks by assuming that the effective masses satisfy 



M{A{BC}) = M{A{DD}) and M{A[BC]) = M{A[DD]) , (7.13) 

where uid = {mB + "ic)/2. We make clear in the following where we are using this 
assumption and where not. 

A large part of our analysis concerns mass splittings between, for example, spin- 1/2 
baryons. In most cases, we have extracted these differences by directly fitting to the ratio 
of the appropriate correlators. Thus, for example, 

^^(^) _ p-(Ms-Mjv)t 

and so the mass difference is obtained from a single fit. This has the advantage of both 
reducing some of the systematic errors (e.g. those arising from excited state contamination) 
and of improving the statistical errors. We find that estimates obtained in this manner are 
consistent with those obtained from individual fits but have errors which are 3 — 5 times 
smaller. 



7.2. Spin- 1/2 baryon mass splittings 

In this subsection we use our results for mass splittings between spin-1/2 baryons to 
extract predictions for the physical mass splittings, and to study the chiral behavior of 
baryon masses. It turns out that the leading chiral prediction, Eq. (|7.9|) , gives a poor 
description of our data. Higher order terms are essential in order to extrapolate reliably 
from our quark masses, and have a significant impact on the final result. The higher order 
terms we use are motivated by quenched chiral perturbation theory — which predicts non- 
analytic terms of 0(mg^^) and analytic terms proportional to m^. The form of these terms, 
and some predictions for their coefficients, have been worked out in Ref. We attempt 



only a qualitative comparison with these predictions — though we provide the results of our 
fits so that others can pursue this further. 

We first consider the octet hyperfine splitting, "E — A" , 

Ma{bb} - Ma[bb] ^ (_8^/3) ^ ^ Mls-M%s ^ ^. Mj^-Mj^ 

niA-rnB niA - ms ^ niA - ms (7.14) 

+ einiB + e[{mA + ms) , 

where Mab is the mass of the pion with fiavor AB, etc. The constants di and can be 
expressed in terms of parameters of the quenched chiral Lagrangian. For reasonable choices 
of these parameters, one expects \d'i\ <^ \di\. There is no useful information concerning ei 
or e'l. 
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0.2 0.4 0.6 



Fig. 9: Test for chiral corrections in My: ~ M/^. The value extrapolated to 
the physical point is shown by the burst symbol at the extreme left. 
In this, and subsequent, figures all masses are in lattice units. 



We fit our results in two ways. First, we assume that the di term is dominant, and 
plot our data versus (M^g — M^^)/ (myi — m^). The outcome is shown in Fig. ^ The data 
should collapse onto a single curve, which should be linear, and our results are reasonably 
consistent with this. What is particularly striking, however, is the size of the slope. This 
is a clear sign that terms of higher order than linear in the quark mass are needed to 
describe our baryon masses — recall that the linear term has been divided out in this fit. 
If we extrapolate linearly to the physical point (m^ = mg, ms = m) we find (the fit gives 
D = -0.36(3), di = -0.19(2) GeV"^) 

Me -Ma = 76(7) MeV. (7.15) 

It is conventional to quote this result in terms of an effective D parameter 

D,, , . Z!2P)^ . ^0.28(3) . (7.16) 
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Here (and in the similar estimates of D^s and F^s below), we use the non-perturbative 
definition of the MS quark mass calculated in MeV, but depart from our usual practice 
by evaluating it at the scale 2 GeV (instead of 1/a). Note that -Deff differs from the D 
appearing in chiral expansions such as Eq. ( [7.14|) , for it absorbs some of the higher order 
terms. 
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Fig. 10: Quadratic fit to Mj^ — M^/ (jn — nis), including baryons composed of 
completely non-degenerate quarks. The value extrapolated to the 
physical point is shown by the burst symbol at the extreme left. 



In our second fit we test to see whether our data can be represented as well by the 
analytic terms. We have found, by trial and error, that Eq. ( |7.14| ) with di = d'l = 
and ei = e'l does a reasonable job; the plot is shown in Fig. |10|. In this plot we have 
included the results for baryons composed of completely non-degenerate quarks, based on 
the discussion following Eqs. ( 7.11|) and (|7.12|) . To emphasize this, we have labeled the 
X-axis by -|- + nis instead of tua + 2mB. The data show definite curvature, so we 
have extrapolated to the physical point using a quadratic fit (corresponding to terms up 
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to in the original expression for baryon masses!). The fit yields 



Me - Ma = 80(8) MeV , D^s 



■^"(^'^ = -0.29(3). (7.17) 



{tUs — to) 



3/2 

These results are consistent with those from the "to^ " fit. 



2.5 



m 



1.5 




0.2 



0.4 

r3 TV/-3 



0.6 



Fig. 11: Test for chiral corrections in M-^ — M^r as described in Eq. ( |7. Jg| ). 

Tie physical point is shown by the burst symbol at the extreme 
left. 



Next we consider the "S — N" splitting, i.e that between two E-like states having 
different strange quark masses. The chiral expansion is 

MAjcc} - Ms^cc} _ ^ - ^ ^, Ml^ - M|^ 

niA-mB uiA-m^B ^ m,A - ms (7.18) 

+ e2TOc + e'2{mA + tob) . 

Chiral perturbation theory suggests that [(iy < Msl- Thus we first plot the data assuming 
d2 is the dominant coefficient (Fig. pT|) . The collapse onto a single curve is reasonable, and 

28 



2.5 



1.5 








0.05 0.1 

(mi + m2) + av(m3) 



0.15 



Fig. 12: 



Quadratic fit to — Mj^ /{m^ — ) for tlie 30 quark combina- 
tions. Tiie result of the two ways of extrapolating to the physical 
value described in the text are shown by the burst symbols. 



a linear fit (F - D = 1.30(5), ^2 = -0.31(3) GeV"^) gives 



Me - Mat = 227(14) MeV, Fes - DefF 



114(7) MeV 

(ms — m) 



1.11(7). 



(7.19) 



We also have investigated various analytic fits, none of which do a good job of collaps- 
ing the data onto a single curve. Our best attempt, shown in Fig. 12, assumes d2 = d2 = 



and 62 = 62/4, i.e. we plot against the average mass of the quarks in the two baryons. We 
again include non-degenerate baryons in this plot. Since the collapse is not good we use 
two different extrapolations to the physical point. The first is a quadratic fit in the average 
mass as shown in Fig. |T^. In the second we first extrapolate linearly to ms = rnc = fn 
and then interpolate to mA = rns- These two methods give almost identical results, as 
shown by the two "bursts" in the plot. The mean of these two points gives 



Me - Mat = 260(20) MeV , F^s - D^s 



130(10) MeV 

(mg — m) 



1.27(10) . (7.20) 
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The difference between the two estimates, Eqs. ( [7.19|) and ( [7.20|) , is indicative of the fact 
that neither of the forms fits all the data well. 




0.2 0.4 0.6 

(Mab - Mlc)/(mB - mc) 

Fig. 13: Test for chiral corrections in M~, — Mjy as described in Eq. ( \7.21\ ). 

The physical point is shown by the burst symbol at the extreme 
left. 



Thirdly, we consider the difference "S — A^" , in which we study the effect of changing 
the mass of the two symmetrized quarks. We expect 

niB -mc TTiB -mc ^ niB - m^c (7.21) 

+ esm^A + e'^imB + mc) . 

It turns out that in quenched chiral perturbation theory we expect d^, = d2 and 63 = 62. 
In this case, there is no expectation that d^ and d'^ should be substantially different in 
magnitude. Nevertheless, for simplicity, we first plot the data assuming (is is the dominant 
coefficient — as shown in Fig. O. The data collapses reasonably onto a single curve, and a 
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Fig. 14: Quadratic fit to Me - MN/{m% - mf„ 

Tiie result of the quadratic extrapolation to the physical value is 
shown by the burst symbol. 



linear fit (F = 0.92(3), rfg = -0.30(3) GeV"^) yields 

M= - Mtv = 338(19) MeV, Feff = = 0.82(5). (7.22) 

(rris — m) 

A quadratic fit to the average quark mass (i.e. assuming 63 = e'g) is slightly better, as 
shown in Fig. ^ The results are nevertheless consistent with those from the first fit: 

M= - Mat = 355(21) MeV, F^s = ^^^^^ = 0.87(5). (7.23) 

(ms — m) 

It is important to note that, once again, higher order corrections in the chiral expansion 
are substantial. 

To further investigate the mass splittings we have considered the combination of dif- 
ferences appearing in the Gell-mann-Okubo formula, 3Ma + M-£ — 2Mn — 2M=; = 0. In 
our notation, this combination is 

GMO(A, B) = 3Ma[bb] + Ma{bb} - ^Mbibb} - 2Mb{aa} ■ (7.24) 
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It is constructed to cancel the 0{mq) terms, 



GMO(A, B) = d^{MlA - 2M1b + MIb) + e^{mA - mef + ■ ■ ■ , (7.25) 

making it an interesting window on higher order terms. We have calculated GMO{A, B) 
using three methods: 

(1) simply taking the differences of the masses; 

(2) using 

GMO(A, B) = A[Ma{bb} - Mb{bb}] - ^Ma{bb} - Ma[bb]] ^ ^ 

(7.26) 

- 2[Mb{^^} - Mb{bb}] , 

where the differences within square brackest are calculated from ratios of correlators; 
and 

(3) using another combination of differences (calculated using ratios) 

GMO(A, B) = [Ma{bb} - Mb{bb}] + ^Ma[bb] - Mb{bb}] 

(7.27) 

- 2[Mb{aa} - Mb{bb}] ■ 



Our results are given in Table ^ They are consistent with zero for all quark mass 
combinations, showing that higher order chiral corrections are not uniformly large. This 
is in qualitative agreement with experiment, where the GMO relation works well: GMO = 
26 MeV, or GMO = 0.011 in lattice units. It is difficult to make a quantitative comparison, 
since we cannot extrapolate our data to the physical point. We note, however, that our 
results for the largest mass splitting, A/B = S/Us, are consistent with the experimental 
value. 



Fit: 


3/2 

niq' 


ml 


Eq. O 


Expt. 


Me - Mn 


227(13*) 


260(20*) 


194(12*) 


253 


Ms - Mn 


338(19*) 


355(21*) 


306(22*) 


375 


Ms-Mj: 


107(9*) 


109(8*) 


112(14*) 


122 


Me - Ma 


76(7*) 


80(8*) 




77 



Table 13. Estimates of mass splittings in the baryon octet, using various fits explained in 
the text. Experimental results are given for comparison. All results are in MeV. 



A summary of mass-splittings in the octet multiplet is given in Table |TB|. Several com- 

3/2 

ments are in order. First, the "mg " fits show evidence for curvature which, if included, 
would likely make the results agree more closely with those from the "mg" fits. Second, we 
reiterate that the data is extremely poorly represented by the first order mass formulae, 
Eq. ( [7.9| ) — this would predict that all the curves in Figs. |9|-|I^ are flat. If, nevertheless. 
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Table 12. Tests of tiie Gell-mann-Okubo mass formula: results for GMO(A, B). 



we fit to the linear terms, we find the results listed in the third column. Finally, we note 
the good agreement of the results from the fit (which are our most reliable) with the 
experimental splittings. 

To further investigate the applicability of quenched chiral perturbation theory we have 
constructed several other quantities. Although most of these are peculiar to the quenched 
approximation, having no counterpart in QCD, they allow us to see how well we understand 
the extrapolations that are needed for all quantities, including those which are physically 
relevant. 

1. We begin with the double difference 

Xe(A, S, C) _ Ma{cc} - Mb{cc} - Mc{AA} + Mc{bb} 



niA — TUB TnA — ITlB 

rriA - ttlb 



— 



(7.28) 

This is interesting for several reasons. First, the expected chiral form is simpler than 
the differences considered above. Second, for C = S or C = A, Xy. reduces to 
M{A{BB]) - M{B{AA}), which for A = s and S = w is rriE. - mj:. Third, is 
predicted to be independent of m,c- Strictly speaking, this is true up to quenched 
artifacts proportional to S. Thus the mc dependence of X-^ is a window onto such 
artifacts. Our data for Xj:, plotted in Fig. Il5|, confirms our expectations. There is no 
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significant dependence on mc — indeed, there is barely any dependence on or ms 
either. The fit yields F + D = 0.50(6) and 4 - 4 = -0.013(16) GeV"^ assuming 
^2 ~ ^3 = 0- Clearly we could just as well fit with the e' terms. If we extrapolate to 
the physical point (assuming M, 



2 

ss 



2M| 



M^), we find 



53f5l MeV 

Ms -Me = 107(9) MeV, {F,s + D,s) = / ^ _, = 0.52(5). (7.29) 



consistent with the other estimates given in Table |T^. It is noteworthy that the ex- 
trapolation required is minimal, unlike that for most of the mass differences considered 
above. 
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Fig. 15: The quantity Xj: determined from Gts to the ratio of correlators. 

The extrapolated value from which M~, — Me is calculated is given 
by the burst symbol at the left. 
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2. We next consider the corresponding quantity for the A-hke baryons, 



XAiA,B,C) ^ Ma[cc] - Mb[cc] - Mc[aa] + Mc[bb] 

TUa — rriB TTlA — TTlB 

= -2F + lOD/3 + d5^^ ^ + e^iruA + tub] 

rriA - tub 



(7.30) 



This always involves a quenched particle with no correspondent in nature, e.g. 
M{U[SS]). We plot Xa in Fig. Again the expectation of no dependence on 
mc is borne out, but this time there are significant higher order terms. The fit yields 
2F - 10L'/3 = 3.4(3) and 4 = 0.26(4) GeV"^ assuming 65 = 0. The data is equally 
well described using a 65 term with ^5 = 0. 
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Fig. 16: The quantity Xa determined from the masses of the individual 
states. The point for physical quark masses is the shown with the 
burst symbol at the left. 



3. In the quenched approximation, we can form two additional double differences, anal- 
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ogous to GMO(A, B), in which the 0{mq) terms cancel 



Y^{A,B,C,D) 
Ya{A,B,C,D) 



M 
M 



A{CC} 



M 



B{CC} 



M 



A{DD} 



+ M 



B{DD} : 



A[CC] - Mb[CC] - Ma[DD] + Mb[DD] ■ 



(7.31) 



Yy, picks out the ^2 ^'2 terms in Eqs. (|7.18|) and ( |7.28| ). The predicted form is 
Yj:{A,B,C,D) {ijiA- mB){mc - mo) 



M'ac - M%c 



d2 + 62 



(7.32) 



We test this in Fig. |T^. What is most noteworthy is that, unhke the GMO re- 
lation, there are significant higher order corrections. We find a reasonable fit if 
d2 ~ —5.4(1.4) GeV~^ and 62 ~ 44(11) GeV~^. The large errors in the fit pa- 
rameters are due to the fact that the two variables M\q — M^c ~ ^ad + ^bd 
{mA — 1^3) (i^c — ^d) are nearly proportional for our set of quark masses. There is 
thus a significant cancellation between the two terms in the fit. Indeed, we can obtain 
as good a fit setting 62 = 0, and replacing it with a higher order term proportional to 
rriA + mB + mc + tod- 
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Fig. 17: Testing the mass dependence of Y-^ (itself determined from fits to 
the ratio of correlators). 
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4. Finally, we consider 1a, which is expected to have the same functional form as Yj2. 



The data, shown in Fig. ^ is again consistent with the expected form, although the 
errors are larger since I^a is calculated from masses instead of mass differences. The 
fit gives -4.3(2.4) GeV"^ for the intercept and 35(20) GeV"^ for the slope. 
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Fig. 18: Tiie quantity Ya determined from the masses of the individual 
states. 



In summary, we have demonstrated that terms beyond linear order in the quark mass 
are necessary to fit the octet baryon mass splittings. The lever arm provided by our 
four "light" quark masses is not, however, strong enough to disentangle the m'^/^ and 
contributions. Indeed, it is important to note that, if we did not have results for the 
baryons containing the S quark, our evidence for higher order terms would have been 
much weaker, and our data would have been reasonably well fit by an 0(m) term alone. 
This would, however, have led to underestimates of all the mass differences, particularly 
mj: - m,N- 

7.3. Nucleon mass 

We now turn to the overall mass scale of the spin-1/2 baryons, which we set using the 
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Fig. 19: Linear fit to and Ma for the three degenerate quark combina- 
tions. The fourth degenerate point SSS is also shown. 

nucleon mass. As for the mass splittings, we expect to need terms of higher order than 
linear in the chiral expansion to describe our data, so we analyze the data with and without 
them. We first fit to the masses of the baryons containing three degenerate quarks. The 
data are shown in Fig. |T^ (along with those for the spin-3/2 baryons). We begin with a 
linear fit to the lightest three baryons, yielding the fit shown in the Figure. We find 

Mn = (0.461(9) + 4.7(2)m) = 1084(27) MeV. 

3 /2 

Next, we fit all four masses including an niq term, yielding Mn = 1070(35) MeV. Finally, 
we fit all four masses including an term, which gives Mn = 1072(31) MeV. In the last 
two cases we have included the same form for the higher order corrections in Mp for the 
determination of the scale and m. Thus the inclusion of curvature systematically reduces 
Mn/Mp by about O.Scr, although the precise functional form of the higher order terms 
is not resolved. For our best estimate we take the mean of the last two estimates, i.e. 
Mn = 1071(35) MeV, corresponding to a'^^MN) = 2062(56) MeV. 

The analysis of the previous paragraph is based on our full sample of lattices. Most 
of our results for mass differences came, by contrast, from our sub-sample of 110 lattices. 
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Thus it is interesting to repeat the extraction of m^r on this sub-sample. We find that 
the results are slightly lower, though consistent within errors. For example, the linear 
"3 point fit" to the lightest degenerate baryons yields M^r = (0.446(9) + 5.04(21)m) = 
1051(26) MeV, about la below the corresponding result from the full sample. Using the 
sub-sample, however, we can attempt a global fit including baryons with non-degenerate 
quarks. We have done this using the S-like baryons composed of Ui quarks, fitting them 
to the linear chiral form of Eq. (|7.9| ), We include all eighteen Ua{UbUc} correlators, using 
the prescription of Eq. (|7.13|) . The fit, shown in Fig. gives 

Mn = 0.452(9*) + 1.91(8*)mA + 3.01(16*)(mB + mc)/2 = 1064(26*) MeV . 

This is in excellent agreement with the 3 point fit. The fact that the fit is reasonable 
underscores the point made above that the need for higher order terms is not apparent 
using the Ui quarks alone. A similar statement holds for the A-like states. 
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Fig. 20: Global Gt to nucleoli data using Eq. ( |7.g|) . The six cases for niA = 
Ui are shown by square symbols, uia — U2 by octagons, and ruA — 
Us by crosses. 
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7.^. Spin-3/2 baryon mass splittings 

The analysis of the masses of the spin-3/2 baryons is more straightforward. As noted 
above, the correlators are labeled {ABC}, and are completely symmetric between the 
three flavors. We can form 20 states with our four masses, and none of these mix with 
each other. 
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Fig. 21: Linear fit to Mdecupiet for the 20 quark combinations. The four 
degenerate cases are shown with an octagon symbol 



As for the spin-1/2 baryons, we first concentrate on mass splittings. The chiral expan- 
sion for the spin-3/2 masses has been worked out for the case of 3 non-degenerate quarks 



2^ , and is (excluding quenched artifacts proportional to 5) 

M{ABC} = Mo + AM + c^iniA + ms + mc) 

+ rff (Mis + + M^a) + d^{M\^ + Mis -f Mlc) 
-f e^iniAmB + msmc + mcmA) + e^(m^ -\- m% + m^) . 



(7.33) 



AM is the decuplet-octet splitting in the chiral limit. In contrast to the spin-1/2 baryons 
the terms of higher order than linear in the quark mass are small. This is shown in Fig 
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in which we give the result of a linear fit to all 20 masses. We have also done linear 
fits to the 10 baryons composed of Ui quarks, and to the 3 lightest baryons composed of 
degenerate quarks. We refer to these three fits as "20-point", "10-point" and "3-point" 
fits, respectively. The results are consistent with one another: 

Mdecupiet = 0.590(16 ) + 3.18(34 )mq 3 point fit (fuU sample) , 

Mdecupiet = 0.578(16*) + 3.49(33*)mg 3 point fit , 

Mdecupiet = 0.578(16*) + 3.50(34*) (m^ + ms + mc)/3 10 point fit , 

Mdeeupiet = 0.580(12*) + 3.42(14*) (m^ + ms + mc)/3 20 point fit , 

(7.34) 

although there is a small systematic difference between the results from the full sample 
and the sub-sample. 

The higher order terms, though small, are nevertheless present. We have investigated 
them using fits to ratios of correlators. We first consider Mq — Ma, which is obtained by 
fitting our data to 

rriA - ms tua - rns (7.35) 

+ 3(ef + e^)(m^ + m^) . 

As before, we can fit our data keeping either the + d^ or the ef' + term. The former 
fit is shown in Fig. |2|, and has parameters = 0.94(12) and d^ + d^ = 0.023(16) GeV"^ 



The slope is small and marginally significant. Extrapolating to the physical point (using 
= 2M|, - M2) yields Mq-Ma= 308(24) MeV. 

We have also done the extrapolation using a more traditional method: we calculate 
M{{AAA}) - M{{BBB}) for B = Ui and A = S or Ui, extrapolate linearly to tub = rn, 
and then interpolate linearly to uia = rng. We refer to this as the "3 point, ratio" method. 

3 /2 

The result is Mn — Ma = 316(18) MeV and is consistent with our niq estimate. The 
result from our subsample, 310(19*) MeV, is consistent. 

We have repeated this analysis for the other two physically interesting mass splittings, 
and collect the results in Table ^ We only quote results from the "3 point, ratio" ex- 
trapolations, since, with such mild deviations from linearity, any reasonable extrapolation 
method using our ratio data gives almost the same result. We also include the results 
of extrapolating using the parameters from the 20-point linear fit. The small difference 
between the two sets of results is due to the higher order terms. Both sets of results are 
significantly smaller than the experimental splittings, and the inclusion of the higher order 
terms makes the disagreement worse. These features are in marked contrast to those we 
saw in the spin-1/2 baryons. 

We can study the higher order terms in more detail by looking at violations of the 
"equal-spacing rule" 

M^AAA} - M^AAB} = M{AAB} - M{aBB} = M{abb} - M{bBB} ■ (7.36) 
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Ratio 


Eq. (I7.34P 




State 


3 point 


20 point 


Expt. 


Ms* - Ma 


107(6*) 


116(7*) 


151 


Ms* - Ma 


215(12*) 


232(14*) 


299 


Ma -Ma 


316(18) 


347(21*) 


436 



Table 14. Estimates of mass splittings in the baryon decuplet. The different methods are 
explained in the text. All results are in MeV. 



This rule holds when keeping up to the linear term in the chiral expansion. Experimentally, 
the three splittings are Mq-Me.* = 137.5 MeV, Me.* -My.* = 148.1 MeV and Me* -Ma = 
151 MeV. The violations of the rule are thus small, and it is interesting to see how well the 
quenched approximation reproduces the magnitude and pattern of these violations. We 
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consider two double differences 



ES1(A, B) = {M^AAA} - M^AAB}) - {M^ABB} - M^bBB}) , 
ES2(A, B) = {M^AAA} - M^BBB}) - ^{M^AAB} - M^ABB}) ■ 

The first becomes {Mq — M^*) — {M-^* — Ma) = —13 MeV when uia = Ws and = m. 
The expectation from quenched chiral perturbation theory is that 

ES1(A, B) = 2df{M\A + - 2M1b) + 2ef (m^ - msf , (7.38) 

i.e. the form of the higher order terms is the same as that appearing in the GMO relation 
(Eq. (|7.25| )). They are shown in Fig. |2^, and are consistent with the expected chiral form 
(df = 0.46(43) GeV"^, ef = -8(7) GeV"^), with large errors.. Extrapolating to the 
physical point, we find ESI = —4(7) MeV, marginally inconsistent with the experimental 
value. 




0.22 0.24 0.26 0.28 0.3 0.32 
(m^ - mB)V(ML+ML-2M^B) 



Fig. 23: Test of violations of the equal spacing rule, Eq. ( \7.3di) . The physical 
point is indicated by the burst symbol, here at the right. 
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The second double difference, ES2(A, S), is predicted to vanish by the general form 
of Eq. ( [7.33|) . The same is true in chiral perturbation theory in full QCD ||2^. Thus it is 
a window onto the quenched artifacts proportional to 5, whose form is predicted to be 

ES2(A, B) = 5/^chiral(A, B) , 

chi.al(A B) . - ™«) (ln(M3,/M|«) + 2 - . 

As explained in section Sec. ^, we have omitted such artifacts from previous expressions, but 
we include them here so as to have a form to fit with. Note that the function chiral(A, B) 
is antisymmetric under A <-> S, as it must be in order to match the antisymmetry of the 
l.h.s.. Note also that it diverges logarithmically when uia — an example of the sickness 
of the quenched approximation in the chiral limit. 



Our results, shown in Fig. are statistically different from zero. We show a linear fit 
of ES2(A, i?)/chiral(A, B) versus + m^, which yields an intercept of 5f^ = 0.29(19) 
and a slope of —0.9(0.7) GeV~^. Interpolating to the physical point, we find ES2 = 
—29(18) MeV. This has the same sign as the experimental value 

Ma - Ma- 3(Ms. - Ms*) = -8 MeV , 

but its magnitude is much larger. It is important to realize, however, that our values 
for ES2 range from —1 to —10 MeV. The large value after extrapolation is due to the 
divergence of chiral(A, B) as m ^ 0. We have attempted analytic extrapolations, using a 
linear combination of (m^ — m^)^ and (m^ — mB)(m^ + m^)^, but this ansatz does not 
fit our data. 

The most important conclusion to be drawn from this analysis is that the decuplet 
splittings are smaller than the experimental values. This discrepancy is only made worse 
if we use other definitions of the strange quark mass, e.g. ms{Mx)- We have also found 
some evidence for artifacts due to quenching in the baryon spectrum. Given the size of 
our errors, however, this result is far from definitive. 

7.5. Ma and Ma - Mn 

We close this section with our results for the overall scale of the spin-3/2 baryons, and 
its relation to that of the spin-1/2 baryons. Our preferred value for Ma comes from the 
3-point linear fit to the full sample (shown in Fig. ^9]), yielding Ma = 1382(36) MeV. Fits 
to the sub-sample, the parameters of which are given in Eq. ( [7.34| ), yield values that are 
consistent with one another, but lie slightly below that for the full sample. For example, the 
20-point fit gives Ma — 1362(36*) MeV. We have also tried 4-point fits to the degenerate 
states, including m'^/^ (or m^) corrections in both A and p. These increase the estimates 
of Ma to 1410(48) (1395(43)), the bulk of the increase coming from the change in scale 
due to the curvature in Mp. 

We have calculated Ma — Mn from the ratios of correlators using the full data sample. 
Using linear extrapolation from the three light mass points we find 318(30) MeV. Including 

an mq^"^ or term in both Ma-M^ and Mp fits we get 365(44) and 347(39) respectively. 
These values are slightly higher than the experimental value 293 MeV. A summary of our 
baryon mass results (without extrapolation to a = 0) is shown in Fig. |2^ along with the 
experimental data. 
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Fig. 24: Test of violations of the equal spacing rule, Eq. ( \7.3f^) . The physical 
point is indicated by the burst symbol (second from left). 



We have also made fits to the negative parity baryon states. The signal in the corre- 
lators is much poorer, falling below the noise by t = 12. For this reason we only present 



the summary shown in Fig. 26 



8. Infinite volume continuum results 



There exist three other high statistics calculations of the spectrum with Wilson 
fermions at /3 = 6.0 on lattices of size 24^ and 0. Their results are given in 

Table 15 along with our best estimates. The data indicate that there are no significant dif- 
ferences between the 24"^ and 32"^ lattices. Thus we do not apply any finite size corrections 
to our data. 
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Fig. 25: A comparison of the baryon masses obtained in this calculation with 
the experimental data indicated by the horizontal lines. The scale 
a is set by Mp and the data for all states are from the sub-sample 
of 110 lattices. 



Our best estimates aX (3 — 6.0 for various mass ratios of interest are 



Mat 
Mp 

Ma 
Mp 
Ma 

Mn 



1.412(35) 
1.800(47) 
1.275(36) 



Expt: 1.22, 
Expt: 1.60, 
Expt: 1.31. 



(8.1) 



These are obtained by extrapolating individual masses linearly to m, and then taking ratios 
within the jack-knife procedure. All extrapolations are done using only the Ui quarks. To 
check for extrapolation errors we have also calculated the ratios for each quark mass and 
then linearly extrapolated these to m. This yields Mu/Mp = 1.42(3) Ma/M^ = 1.79(4), 
and Ma/Mm = 1.25(3) consistent with the first method. 
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Fig. 26: A comparison of the negative parity baryon masses obtained in this 
calculation with the experimental data indicated by the horizontal 
lines. Shaded bands show the experimental uncertainty. The scale 
a is set by Mp and the data for all states are from the sub-sample 
of 110 lattices. 



Our results for M^/Mp and M/^/Mp are larger than the experimental values. It has 
been stressed by the GFll collaboration, however, that these ratios have a significant 
dependence on a. In an attempt to check this we combine our results with those from 
the GFll Collaboration 0] at /3 = 5.7, 5.93, and 6.17. The two calculations have similar 
statistics, while the physical volume of our lattices is larger. The GFll collaboration 
have made two separate extrapolations to the continuum limit using different gaussian 
smeared sinks. In their notation, "012" is a combination of results from three sinks with 
smearing radii 0, 1, 2, while "4" refers to the use of a single larger smearing radius of size 
4. In the final analysis they prefer to use the "012" result as it has smaller statistical 
errors. We update these two linear fits (using the same 3 data points used in [|2| and 
ours at /9 = 6.0), and the results are shown in Fig. |2^. After including our point, the 
extrapolated values for the "012" data (dotted lines) change from M^/Mp = 1.28(7) to 
1.30(6) and from Ma/Mp = 1.61(8) to 1.62(7). The x^/dof for the new fits are 2.1 and 0.85 
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0.298(2) 


0.428(4) 


0.647(6) 
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0.296(1) 
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iUUU 


U.zyD4(^4J 


M /I OQ /'0^ 

(jAZdyZ ) 


\j.o4Z[6 ) 




0.1558 


24^ X 32 


78 


0.234(3) 


0.397(5) 


0.574(8) 


0.686(25) 


0.1558 


32^ X 64 


170 


0.234(1) 


0.387(3) 


0.571(4) 


0.664(8) 


0.1563 


243 X 32 


78 


0.184(3) 


0.378(10) 


0.522(14) 


0.636(45) 


0.1563 


243 X 54 


200 


0.184(2) 


0.343(14) 


0.530(16) 


0.629(19) 


0.1563 


32^ X 64 


170 


0.185(1) 


0.361(5) 


0.525(7) 


0.636(13) 



Table 15. Comparison ofhadron masses obtained using Wilson fermions at (3 = 6.0. Results 
from 32"^ x 64 lattices are from the present work, while 24^ x 32 data is from 
Ref pW, 24^ X 54 data is from Ref pW, and 24^ x 64 data is from Ref jO/. 



respectively. The analogous numbers for the sink "4" data (solid line) are 1.33(9) —>■ 1.38(7) 
and 1.68(10) 1.73(10) with / dof for the new fits equal to 1.2 and 0.86 respectively. 

As is evident from Fig. the main difference between the two fits comes from the 
difference in the "012" and "4" data at /3 = 5.7. On the basis of x^/dof^ we find that 
combining our results with the sink "4" GFll data is preferred, in which case there is 
very little a dependence. If we neglect the point at strongest coupling, (3 = 5.7, then the 
remaining three points again show no clear a dependence for both M^/Mp and M^/Mp, 
and give very similar values for the fit parameters. 

The ambiguity in the extrapolation makes it clear that data at more values of {3 are 
needed in order to reliably extrapolate to the continuum limit. Our preferred estimates 
are Mn/Mp = 1.38(7) and Ma/M^ = 1.73(10) from fits shown by a sohd line in Fig. 1^. 
This suggests that the quenched approximation is good to only ~ 10 — 15%. 



9. Conclusions 

We have presented a detailed analysis of the hadron spectrum in quenched QCD 
at /3 = 6.0 with Wilson fermions, focusing on states composed of light quarks. Our 
small statistical errors, and our use of several moderately light quarks, have allowed us 
to improve the extrapolations to the chiral limit. This is particularly true of the mass 
splittings amongst the octet baryons. Here we find substantial contributions from terms 
of higher order than linear in the quark mass. Motivated by quenched chiral perturbation 
theory, we have found good variables with which to extrapolate to the physical quark 
masses. Our results show that the splittings are larger than previously thought, and are 
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Fig. 27: Linear extrapolation to the continuum limit of the ratios Mj^/Mp 
and M^/Mp. Our data is shown with symbol octagon, and the rest 
of the points are from the GFll Collaboration The ^'012" sink 
points are labeled by squares and the sink "4" points by crosses. 
The extrapolated values are shown by symbol fancy cross and are 
shifted from Mpa = for clarity. At jJ = 5.7 (Mpa ~ 0.56) we 
have shown GFll collaboration's data for both 16^ and 24"^ lattices 
even though they used only the 16"^ data in the fit and shifted their 
continuum result by the difference to account for finite size effects. 



comparable with their experimental values. These results emphasize the importance of 
calculating masses of baryons composed of several combinations of non-degenerate quarks. 

The extrapolations required for the tt, p, N and A are less sensitive to higher order 
terms. There is a clear curvature in the p and nucleon channels, and it can be accommo- 

3 /2 

dated either by including a term of 0{mq ) (which would result from chiral loops) or a 
term of 0{m'^). The effect on the extrapolated values for the Mp and Mn is, however, 
small — roughly comparable to the statistical errors. Higher order terms are also small for 
the decuplet baryons. 
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It is not surprising that higher order terms are needed when considering quarks with 
masses ranging up to and beyond that of the physical strange quark. In previous cal- 
culations of baryons composed of degenerate quarks, these terms were small and often 
neglected. What is striking is how this is not true for many of the mass differences involv- 
ing baryons composed of non-degenerate quarks. 

One caveat concerning the results for mass splittings is the fact that there are substan- 
tial systematic errors in the extraction of m^. Different methods lead to results differing 
by up to ~ 20%. This is presumably an error due to quenching, although some part of it 
could be due to discretization. Our favored choice for is that determined by matching 
to the ratio M^/Mp. This gives the largest estimate for m^. 

How does the spectrum at /3 = 6 compare with experiment? We find that the ratio 
Mjs[/mp = 1.41(4) is too high, while M/\/Mjs[ = 1.27(4) is consistent with experiments. 
Using the larger estimate of mg, the octet baryons splittings agree with experiment, while 
those in the decuplet are too small by 30%. 

Of course it is quite possible that there is substantial variation in some of these 
ratios as we extrapolate to the continuum limit. Combining our data with that of the 
GFll collaboration, however, we conclude that there remains considerable uncertainty in 
this extrapolation. Our preferred extrapolation gives M^/Mp = 1.38(7) and M/\/Mp = 
1.73(10), but the systematic errors exceed those from statistics. Thus, in our view, it 
remains an open question how well the quenched approximation represents full QCD when 
extrapolated to the continuum limit. The errors could well be as large as ~ 10 — 15%. 
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